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Abstract A search for supersymmetric particles in fi- 
nal states with zero, one, and two leptons, with and 
without jets identified as originating from 6-quarks, in 
4.7 fb~-^ of v^s = 7 TeV pp collisions produced by the 
Large Hadron Collider and recorded by the ATLAS de- 
tector is presented. The search uses a set of variables 
carrying information on the event kinematics transverse 
and parallel to the beam line that are sensitive to sev- 
eral topologies expected in supersymmetry Mutually 
exclusive final states are defined, allowing a combina- 
tion of all channels to increase the search sensitivity. 
No deviation from the Standard Model expectation is 
observed. Upper limits at 95 % confidence level on fidu- 
cial cross-sections for the production of new particles 
are extracted. Results are interpreted in the context 
of the constrained minimal supersymmetric extension 
to the Standard Model and in supersymmetry-inspired 
models with diverse, high-multiplicity final states. 

Keywords SUSY • Razor 



1 Introduction 

One of the most promising extensions of the Standard 
Model, supersymmetry (SUSY) [l}|9], has been the tar- 
get of a large number of searches at the LHC. Prompted 
by the large predicted production cross-section of coloured 
SUSY particles (sparticles), ATLAS and CMS have per- 
formed inclusive searches for strongly produced squarks 
and gluinos, the superpartners of quarks and gluons flOL 



discriminate the signal processes from the background 
typically include requirements on the missing trans- 
verse momentum (E!^^^^) and the scalar sum of trans- 
verse momenta of all selected physics objects (Ht) plus 



the scalar 



This paper presents a search for strongly produced 
sparticles that makes use of a variety of final states in- 
cluding high transverse momentum jets and zero, one, 
or two leptons (electrons or muons). The events are also 
separated according to the presence of a jet identified 
as originating from a 6-quark (6-tagged jets). Several 
mutually exclusive search channels are defined, facil- 
itating a simultaneous search in all of the typical fi- 
nal states and increasing the search sensitivity. The 
search employs a set of observables, called the "razor 
variables" 



17 . Assuming R-parity conservation 18 -22 , these spar 



tides are produced in pairs and decay into energetic 
jets, possibly leptons, and the lightest SUSY particle 
(LSP, typically the lightest neutralino Xi), which es- 
capes detection and results in missing transverse mo- 
mentum. For these searches, the selections adopted to 
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which make use of both longitudinal 
and transverse event information. These variables were 
first employed in SUSY searches by CMS (24||25]. Their 
use makes this search sensitive to different regions of 
kinematic phase space relative to the other iJip^'^-based 
searches. Thus, these search results complement those 
already performed by ATLAS. 

This paper is organised as follows. The main fea- 
tures of the ATLAS detector are presented in Sect. [2] 
Section [3] introduces the razor variables. Section [H de- 
scribes the data sample, basic event selection, and the 
Monte Carlo simulation used to model the data. Sec- 
tion [5] defines the basic physics objects and event-level 
variables that are used through the analysis. The search 
technique is described in Sect.|6j and the background es- 
timation is presented in Sect. [7] The performance of the 
search and interpretation of the results are presented in 
Sect. [8] Finally, Sect. |9] includes a summary of the anal- 
ysis and of its findings. 
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2 ATLAS detector 



The ATLAS detector comprises an inner tracking de- 



tector, a calorimeter, and a muon system 26 . The in- 



ner detector includes a silicon pixel detector, a silicon 
microstrip detector, and a transition radiation tracker. 
It is immersed in a 2 T axial field and precisely mea- 
sures the tracks of charged particles in the pseudora- 
pidity regiorj^ \ri\ < 2.5. The calorimeter covers the re- 
gion \r]\ < 4.9 and is divided into electromagnetic and 
hadronic compartments. The electromagnetic calorime- 
try in the central (|?7| < 3.2) region is provided by liq- 
uid argon sampling calorimeters with lead absorbers. 
In the barrel region (|?7| < 1.4), the hadronic calorime- 
try is provided by scintillator tiles with steel absorbers, 
and the more forward (1.4 < \t]\ < 3.2) region is cov- 
ered by a liquid argon and copper sampling hadronic 
calorimeter. The forward calorimetry (|?7| > 3.2) uses 
liquid argon and copper or tungsten absorbers. The 
muon spectrometer covers |ry| < 2.7 and includes a sys- 
tem of air-core toroidal magnets. A variety of technolo- 
gies are used to provide precision muon tracking and 
identification for jryj < 2.7 and rapid response for trig- 
gering for |ry| < 2.4. 

ATLAS uses a three-tier trigger system to select 
events. The first-level (LI) trigger is hardware-based 
and only uses coarse calorimeter information and muon 
system information. The calorimeter information avail- 
able at the lowest level includes basic objects with rough 
calibration and simple identification of electromagnetic 
objects (electrons and photons) as distinct from hadronic 
objects (jets). The second-level (L2) trigger and event- 
filter (EF) compose the software-based high-level trig- 
ger (HLT), in which full event reconstruction is run, 
similar to that used offline, in order to accurately iden- 
tify and measure objects. The L2 only examines rj/cf) 
regions that triggered the LI. The EF fully reconstructs 
events that pass L2. 



3 Razor variable definitions 

Searches for sparticles in R-parity-conserving scenarios 
generally make the assumption that the sparticles are 
pair-produced and decay subsequently to an LSP that is 
invisible in the detector. The heavy sparticles produced 

^ ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the 
2/-axis points upward. Cylindrical coordinates (r, (f>) are used 
in the transverse plane, (f> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as rj = — lntan(0/2). 



are either the same type of particle (pair-production) or 
are at the same mass scale. Thus, the production and 
visible energy in the decays are fairly symmetric. Most 
analyses make use of the transverse balance of typical 
pp collision events, or exploit the event symmetry in 
the transverse plane. The razor variables attempt to 
also include longitudinal information about the event 
by making several assumptions motivated by the kine- 
matics of the models of interest. 

In the rest frame of each heavy sparticle, called the 
i?-framc, the sparticle decays are symmetric. In an at- 
tempt to reconstruct the primary produced sparticle 
pair, the razor calculation clusters all final-state par- 
ticles into a pair of objects with four-momenta called 
"mega-jets." Each of these mega-jets is associated with 
one of the two SUSY decay chains and represents the 
visible energy-momentum of that produced sparticles. 
All possible combinations of the four- vectors of the vis- 
ibly reconstructed/selected objects (signal jets and lep- 
tons) are considered when constructing the two mega- 
jets. The pair of mega-jets, ji and j2, that minimises 
the sum of the squared masses of the four- vectors is se- 
lected. Following the prescription in Ref. [23j and for 
consistency with Ref. 
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all jets and the mega-jets 
are forced to be massless by setting their energy equal 
to the magnitude of their three momenta. Studies in- 
dicate that neither this choice nor the mega-jet selec- 
tion, based on minimizing the mega-jet mass squared, 
have a significant impact on the reach of the razor-based 
search. 

In the i?-frames, each heavy sparticle should be nearly 
at rest with some mass mnoavy The sparticle decay 
may then be approximated as a two-body decay to 
some visible object (a mega-jet) and the invisible, stable 
LSP. The final visible decay products (i.e. the final-state 
quarks and gluons, or the observable jets and leptons) 
have masses far below the SUSY mass scale and can 
therefore be approximated as being massless. Then the 
energy of each mega-jet in the i?- frame, Ei and E2, 
becomes: 



El — E2 — 



2 2 

"^Hcavy ^ "^LSP 
2 X TTlHcavy 



(1) 



where rriLSp is the mass of the LSP. This leads to a 
characteristic mass, M/j, in the R frame of Mj^ — 2 x 
El = 2 X E2, which for mncavy ^ wlsp is identical to 
mHcavy Therefore, in events where heavy particles are 
pair-produced, Mr, which is a measure of the scale of 
the heaviest particles produced, should form a bump. 
In tt or WW events, for example, the characteristic 
mass Mil ~ ^top or mw- The product of Mjf and the 
Lorentz factor for the boost from the lab to i?-frame. 
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— 7-R ^ -^-^-Rj is useful for characterisation of the 4 Data and Monte Carlo samples 
sparticle mass scale, in part because of its close rela- 
tion to mHoavy, and in part because Standard Model 
backgrounds tend to have small values of M^. When 
expressed in terms of the mega-jet quantities in the lab 
frame, the expression is given by: 



The data included in this analysis were collected be- 
tween March and October 2011. After basic trigger and 
data quality requirements, the full dataset corresponds 
to 4.7 ±0.2 fb" 
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■J2,z) 



(2) 



where ji^E and ji^z are the energy and longitudinal mo- 
mentum, respectively, of mega-jet i. The transverse in- 
formation of the system is taken into account by con- 
structing a transverse mass for the mega-jets, assuming 
half of the E!f^^^ is associated with each jet: 



X |E™'""| 



X E™'"" 



X UJl.T 
(jl,T + j2,T) 



|j2,T|) 

1/2 



(3) 



where E™'^^ is the two-dimensional vector of the E™'^^ 



in the transverse plane. When an event contains "fake" 
^miss fj-Qjj^ detector defect or mismeasurement, the 
system will tend to have back-to-back mega-jets. In 
such cases, the vector sum of the two mega-jet momenta 
will be small. If, on the other hand, there is real E^^^^, 
the mega-jets may not be back-to-back and may even 
point in the same direction. In these cases, the vector 
sum, and thus , will have a large value. M-p is an- 
other measure of the scale of the event that only uses 
transverse quantities in contrast to longitudinal quan- 
tities in M'^f . 

Finally a razor variable is defined to discriminate 
between signal and background: 



R 



(4) 



This variable takes low values for multijet-like events 
and tends to be uniformly distributed between and 1 
for sparticle decay-like events, providing good discrimi- 
nation against backgrounds without genuine E^^^^^. The 
impact of some important experimental uncertainties, 
like the jet energy scale uncertainty, are reduced in this 
ratio. In an analysis based on the razor variables, a cut 
on R can be used to eliminate these backgrounds be- 
fore a SUSY search is made in the distribution of the 
variable Af^. 



Events in the zero-lepton channels are selected using 
a trigger that requires a jet with transverse momentum 
PT > 100 GeV at LI. In the event filter, Ht > 400 GeV 
is required, where i/x is calculated through a scalar sum 
of the pt of all calorimeter objects with > 30 GeV 
and 1 77 1 < 3.2. With the exception of a cross-check of 
the niultijet background estimate, which uses prescaled 
single-jet triggers, this trigger requirement is fully effi- 
cient for the offline selection used in the analysis. 

The one- and two-lepton channels make use of the 
lowest-px single-lepton triggers available for the entire 
running period. The muon triggers require a muon with 
Pt > 18 GeV, and the electron triggers require an 
electron with px > 22 GeV. Offline, the leading lep- 
ton in the event is required to have px > 20 GeV 
(px > 25 GeV) if it is a muon (electron) , in order to en- 
sure that the triggers are fully efficient with respect to 
the offline event selection. For the two-lepton analysis, 
where there are overlaps in the triggers, the electron 
trigger takes priority over the muon trigger. 

Offline, an event is required to have at least one 
vertex with at least five tracks associated to it, each 
with p^'^'^^ > 400 MeV. This requirement reduces cos- 
mic ray and beam-related backgrounds. The primary 
vertex is defined as the one with the largest J^iPr'^'^^)'^ 
of the associated tracks. Events that suffer from spo- 
radic calorimeter noise bursts or data integrity errors 
are also rejected. 

Monte Carlo (MC) simulated events were used to 
develop the analysis and assist in estimations of back- 
ground rates. All MC samples are processed through 
atlas's full detector simulation j2^ based on Geant4 
which was run with four different configurations cor- 
responding to detector conditions of four distinct op- 
erating periods of 2011. The fractions of MC simula- 
tion events in these four periods match the fractions 
of data in each period. During the data collection, the 
average number of proton-proton collisions per bunch 
crossing in addition to the one of interest ("event pile- 
up" or simply "pile- up" ) increased from approximately 
two to twelve. To mimic the effect of pile-up, additional 
inelastic proton-proton collisions are generated using 
Pythia 



31 



30 



and overlaid on top of every MC event. 
Within each period, the profile of the average num- 
ber of events per bunch crossing {{fi)) is re- weighted 
to match the data in that period. The same trigger se- 
lection is applied to the MC simulation events, which 
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are then passed through the same analysis code as the 
data. Reconstruction and trigger efficiency scale factors 
are applied to the MC simulation in order to take into 
account small discrepancies between the data and the 
MC simulation. 

Table [1] lists the major backgrounds along with the 
chosen estimation method (described in Sect. [7| and 
the primary and alternative MC generators used in this 
analysis. In all cases, MC@NLO and Alpgen are in- 
terfaced to Herwig and Jimmy for the parton shower, 
hadronisation, and underlying event modelling. The mul- 
tijet background is normalised to the leading order gen- 
erator cross-section predicted by Pythia. The tt pro- 
duction cross-section of 166.8 pb is calculated at ap- 
proximate NNLO in QCD using Hathor 
MSTW2008 NNLO PDF sets 



This model is generated using Herwig H — h 49 , with 



32 with the 



33 . The calculation is 



cross-checked with an NLO-I-NNLL calculation '34] im- 
plemented in Top-|— I- 35 . The single-top production 



cross-sections are calculated separately for s-channel. 



t-channel, and Wt production at NNLO 36 -38 . 

The W and Z (including Drell-Yan with mu 
40 GeV) production cross-sections of 10.46 nb and 0.964 nb 



> 



are calculated at NNLO using FEWZ 39 . For the pro 



duction of vector bosons in association with heavy flavour, 



in accordance with ATLAS measurements 40 , the pro- 
duction cross-section ior W + b and W + cc are scaled 
by 1.63, and the cross-section for W + c is scaled by 1.11 
compared to the NLO cross-section |41 . Additional un- 
certainties on the production of W and Z bosons in 
association with heavy flavour of 45 % for -I- & and 
W + cc, 32 % for W + c, and 55 % for Z + bb are in- 
cluded. Alpgen describes the jet multiplicity and in- 
clusive M'j^ distributions well, but it does not correctly 
model the vector boson distribution. Therefore, the 
boson pt in the Alpgen samples is re-weighted ac- 
cording to the distribution produced Sherpa. Half of 
the difference between the weight and unity is applied 
as a systematic uncertainty on the re-weighting proce- 
dure. Further systematic uncertainties on the shapes of 
Alpgen samples are derived by systematically varying 
the generator parameters, including matching and fac- 
torisation scales. Diboson production cross-sections of 
44.92 pb, 17.97 pb, and 9.23 pb for WW, WZ, and 
ZZ (including off-shell production with m« > 12 GeV) 
are calculated at NLO using MCFM 



the gluino and LSP masses being the only free parame- 
ters. The top quarks are required to be on-shell, limiting 
the mass splitting between the gluino and the LSP to 
greater than 2 x mtop- 

The second considers gluino pair-production, with 
the gluino decaying to two quarks and a chargino via 
an off-shell squark. The chargino then decays to a. W 
boson and the LSP. The free parameters of this model 
are the masses of the gluino, chargino, and LSP. For 
convenience, two two-dimensional planes are generated: 
one with the chargino mass exactly between the masses 
of the gluino and the LSP and one with the mass of the 
LSP fixed to 60 GeV. Because initial-state radiation can 
be important for the acceptance of these models when 
the mass splitting between the gluino and LSP is small, 
this model is generated using MadGraph 50| with at 
most one additional jet in the matrix element. Pythia 
is used for the parton shower and hadronisation. Sys- 
tematic uncertainties on matrix element matching and 
initial-state radiation modelling are included, leading to 
20 % uncertainties for small mass splittings and small 
gluino masses, but no uncertainty for mass splittings 
above 200 GeV and masses above 400 GeV. 

Additionally, the results are interpreted in terms of 
SUSY signal models based on the constrained minimal 



supersymmetric model (CMSSM or MSUGRA) (18 -22 



42 . In order to 



avoid low-mass resonances, all dilepton events are re- 
quired to have the invariant mass mu > 20 GeV. These 
cross-sections provide the starting normalisations for all 
background processes. 

Two SUSY-inspired simplified models are used for 
the interpretation of the results from this search. The 
first considers gluino pair-production, with the gluino 
decaying to a tt pair and the LSP via an off-shell stop. 



The parameters of this model are the high-energy-scale 
universal scalar mass, mg, the universal gaugino mass, 
mi the ratio of the vacuum expectation values of the 
two Higgs fields, tan(/3), the tri- linear coupling strength, 
Aq, and the sign of the Higgsino mass parameter, [i. 
Samples are generated in a two-dimensional grid of the 
mo-mi/2 parameters where tan(/3) = 10 and = 
are fixed and /i is set positive. This MC data grid is 
generated using HerwigH — h 49 , with a more dense 



population of points at low mass. IsaSUSY 51 is used 



to run the high-energy-scale parameters down to the 
weak-scale. 

Signal cross-sections are calculated to next-to-leading 
order in the strong coupling constant, adding the re- 
summation of soft gluon emission at next-to-leading- 

The nom- 



logarithmic accuracy (NLO-I-NLL) [52 56 



inal cross-section and the uncertainty are taken from 
an envelope of cross-section predictions using different 
PDF sets and factorisation and renormalisation scales. 



as described in Ref. 57 . For each of these signal models. 



the luminosity systematic uncertainty of 3.9% 27 28 



and statistical uncertainty, typically of order 10%, is 
included. 
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Table 1 Background estimation methods, primary and alternative MC event generators, and normalisation uncertainties 
for each of the major backgrounds. The backgrounds are constrained using various Control Regions (CRs) that are enriched 
in certain samples (see Sect.|6|. The tt background estimate includes small contributions from ttW and ttZ, generated with 
MadGraph. The diboson WW background estimate also includes W^W^jj generated with MadGraph. The last column of 
the table indicates the uncertainty on the normalization in the simultaneous fit used to test signal hypotheses. "None" indicates 
that the normalization is fully constrained in the fit. The grouping indicates the samples that are collected together in the fit. 
Within a group, the relative normalizations are fixed. 



Background 


0-lepton 


1-lepton 


2-lepton 


Generator 




Alternate 




Normalisation Uncertainty 


Multijets 


MJ CRs 


Matrix method 


Matrix method 


Pythia |311 


Alpgen 






None 


W ^ 


W CRs 


W CRs 


Matrix method 


Alpgen 


43 










None (grouped with Z) 


Z -j. a 


Z CRs 


Z CRs 


Z CRs 


Alpgen 


43 










None (grouped with W) 


Drell-Yan 


Z CRs 


Z CRs 


Z CRs 


Alpgen 


43 










None (grouped with W) 


Z ^ w 


Z CRs 


Matrix method 


Matrix method 


Alpgen 


43 










None (grouped with W jZ) 


tf(had) 


tt CRs 


tt CRs 


tt CRs 


MC@NL 


U 


44] 








None 


^£(leptonic) 


tt CRs 


tt CRs 


tt CRs 


Alpgen 


43 




MC@NLO 


44^47] 


None 


Single top 


tt CRs 


tt CRs 


tt CRs 


MCONL 


U 


44] 








None (grouped with tt) 


WW diboson 


MC 


MC 


MC 


Herwig 


48 




Alpgen 


43 




NLO ± 30 % 


Other diboson 


Z CRs 


Z CRs 


Z CRs 


Herwig 


48 




Alpgen 


43 




None (grouped with WjZ) 



5 Physics object identification and selection 

Events are categorised into six exclusive samples de- 
fined by the presence of zero, one, or two leptons, with 
or without ^-tagged jets. The particle candidate selec- 
tions that define these samples are referred to as the 
"baseline" object selection. Since a particle may simul- 
taneously satisfy multiple particle hypotheses (e.g. elec- 
tron and jet), an overlap removal procedure (described 
below) assigns a unique interpretation to each candi- 
date. The selections are then further refined to define 
"signal" candidates with higher signal-to-background 
ratio in each sample. 

Baseline electrons are required to have E't > 10 
GeV, be within the fiducial acceptance of the inner de- 
tector (I77I < 2.47), and pass a version of the "medium" 
selection criteria (58| updated for 2011 running condi- 
tions, which requires hadronic calorimeter energy depo- 
sition and a calorimetric shower shape consistent with 
an electron and a match to a good quality inner detec- 
tor track. Signal electrons are required to be isolated 
from other objects and satisfy "tight" selections. The 
tight selection applies stricter track quality and match- 
ing than medium and ensures the number of hits in 
the transition radation tracker is consistent with the 
electron hypothesis. The isolation requirement is that 
the sum of the of all charged particle tracks associ- 
ated with the primary vertex within AR = 0.2, where 
AR = a/ (Z\r?)^ + (Z\</))^, of the electron is less than 
10% of the electron i^x- In the leptonic channels, if 
the leading lepton in a data event is an electron, it is 
additionally required to match to an EF trigger elec- 
tron. MC simulation events are re- weighted to compen- 
sate for mis-modelling of the single-lepton trigger effi- 
ciency. The energy of electrons in simulated events is 
also smeared prior to object selection in order to repro- 



duce the resolution in Z and J jt]} data. Finally, in or- 
der to account for differences in electron reconstruction 
efficiency, r/- and i?T-dependent scale factors, derived 
from Z, W ^ and J/V" events in the data, are applied to 
each simulated electron satisfying overlap removal and 
selection requirements. 

Baseline muons are reconstructed as either a com- 
bined track in the muon spectrometer and inner de- 
tector, or as an inner detector track matching with a 
muon spectrometer segment [59| . Tracks are required 
to have good quality, and the muon is required to have 
Pt > 10 GeV and |?7| < 2.4. Signal muons are re- 
quired to be isolated by ensuring that the sum of the 

of all charged particle tracks associated with the pri- 
mary vertex within AR — 0.2 of the muon is less than 
1.8 GeV. Matching to EF trigger muons in data, MC 
event trigger re- weighting, muon momentum smearing, 
and MC/data efficiency scaling are performed in a sim- 
ilar way for muons as electrons (described above) [60]- 



62 



Calorimeter jets are reconstructed from topologi- 
cal clusters of energy deposited in the calorimeter cali- 



brated at the electromagnetic (EM) scale 63 using the 
anti-ftt jet algorithm [64[|65| with a four-momentum re- 
combination scheme and a distance parameter of 0.4. 
Jets reconstructed with an EM-scale > 7 GeV are 
calibrated to the hadronic scale (particle level) using 
and 77-dependent factors, derived from simulation and 
validated with test beam and collision data 66 . In or- 



der to remove specific non-collision backgrounds, events 
are rejected if they contain a reconstructed jet that does 
not pass several quality and selection criteria 66 . Sig- 



nal jets are selected if they lie within |?7| < 2.5 with a jet 
vertex fraction (JVF) of at least 75%, where the JVF 
is the fraction of summed px of the tracks associated 
with the jet that is carried by tracks consistent with 
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the primary vertex of the event, thus associating the 
jet with the pp coUision of interest. Jets are tagged as 
heavy flavour using the combined neural network "jet 



fitter" algorithm 67 with the 60 % efficiency working 



point. Scale factors for heavy flavour jets are used in 
MC simulation in order to reproduce the expected 6-jet 
identification performance in data. 

In order to ensure that objects are not double counted, 
overlaps between objects are removed using a hierachi- 
cal procedure. If any two baseline electrons lie within 
a distance of AR = 0.1 of one another, the electron 
with the lower calorimeter Et is discarded. Next, jets 
passing basic selections are required to be at least 0.2 
units away from all surviving baseline electrons in rj- 
(j). Electrons are then required to be at least 0.4 units 
away from surviving jets. Finally, in order to mitigate 
the effect of jets which have deposited significant energy 
in the muon spectrometer on mass measurements and 
reduce the number of events with badly measured miss- 
ing transverse momentum, muons with px > 250 GeV 
within AR = 0.2 of a jet with pT > 500 GeV are re- 
moved. A negligible number of events in the data are 
removed by this cut. 

Following these overlap removal procedures, the miss- 
ing transverse momentum and razor variables are calcu- 
lated. The determination of the missing transverse mo- 
mentum uses all baseline electrons with Et > 20 GeV, 
all baseline muons, all calibrated jets with pT > 20 GeV, 
and EM scale topological calorimeter clusters not be- 
longing to any object. Note that in the MC simulation, 
objects enter this calculation after the energy or 
smearing described above. 

In counting leptons for event classification, baseline 
electrons and muons are then required to be at least 
0.4 units away from all good jets in ry-^. If an electron 
and muon are separated by ARconc < 0.1, neither is 
counted. 

In order to remove events with large missing trans- 
verse momentum due to cosmic rays, events are vetoed 
if they contain a muon in which the transverse and lon- 
gitudinal impact track parameters are greater than 0.2 
mm and 0.1 mm with respect to the primary vertex, re- 
spectively. Also vetoed are events with badly measured, 
non-isolated muons. These muons with large momen- 
tum uncertainties are rare in both the signal and back- 
ground events and can have significant impact on the 
E!^^^^ and razor variables. 

During a portion of the run period, a hardware fail- 
ure resulted in a region of the calorimeter not being 
read out. For data collected during this period, and for 
a corresponding fraction of the MC samples, events are 



rejected if they fail the "smart LAr hole veto" 13 . This 



ensures that if an event contains one or more jets point- 



ing to the dead region and those jets may contribute 
substantially to the missing transverse momentum in 
the event, the event is discarded. 

Signal regions are defined after all overlap removal 
is complete. Events with no baseline leptons and events 
with the highest-px lepton below the leading lepton re- 
quirement (25 GeV for electrons, 20 GeV for muons) 
are accepted into the zero-lepton regions. Events with 
one leading lepton satisfying all requirements, includ- 
ing that on leading lepton px (above), and no other 
baseline leptons with px > 10 GeV are accepted into 
the one-lepton regions. Events with exactly one addi- 
tional signal lepton above 10 GeV and no other baseline 
leptons are accepted into the two-lepton regions. 



6 Search technique 

After sorting events into the six samples described in 
the previous section, each sample is further divided in 
the R-M'j^ plane into control regions (CR), which are 
choosen so that they are dominated by a specific back- 
ground, and signal regions (SR). Additionally, valida- 
tion regions (VR) are constructed, which do not con- 
strain the background but are used to evaluate the 
agreement between data and MC simulation. Table [2] 
lists these regions, which are also visualised in the R- 
M'^ plane in Fig. [l] These regions are binned in ei- 
ther R or and then simultaneously fit to MC es- 
timates for background and signal rates with corre- 
lations from sample to sample and region to region 
taken into account. The hadronic (had.) and one-lepton 
signal regions are divided into events with and with- 
out 6-tagged jets ("b-tag" and "b-veto," respectively). 
The two-lepton events are divided into regions with 
opposite-sign (OS) and same-sign (SS) leptons and re- 
gions with opposite-flavour (OF) and same-flavour (SF) 
leptons. While some background components are suffi- 
ciently constrained by the CRs to be left free in the fit, 
others are constrained to estimates derived from other 
techniques or MC simulation. Table [T] summarises the 
backgrounds, the estimation technique, the source of 
the estimate, and the normalization uncertainty used 
in the fit. Finally, systematic uncertainties on all back- 
grounds are included as nuisance parameters. The re- 
sult is a maximum likelihood fit that encapsulates all 
knowledge about the background and signal consistently 
across all channels. 

When evaluating a signal sample for exclusion, any 
signal contamination in the control regions is taken into 
account for each signal point, as the control region fits 
are performed for each signal hypothesis. Separately, 
each signal region (one at a time) , along with all control 



Multi-channel search for squarks and gluinos with ATLAS 



7 



regions, is also fit under the background-only hypoth- 
esis. This fit is used to characterise agreement in each 
signal region with the background-only hypothesis and 
to extract fiducial cross-section limits and upper limits 
on the production of events from new physics (-/Vbsm)- 
The fit considers several independent background 
components: 

— tt and single top. A total of five top control regions 
are defined in the one- and two-lepton channels. The 
normalisation of this component is allowed to vary 
freely in the fit. 

— Bosons, except diboson WW . The inclusion of the 
WZ and ZZ diboson samples is motivated by the 
dominance of leptonic Z decays in the two-lepton 
signal regions, the dominance oi Z ^ vv in the 
zero-lepton signal regions, and the dominance of 
WZ — ii/qq in the one-lepton signal regions. In all 
of these cases, the experimental uncertainties affect 
the samples in the same way as they do IV-l-jets or 
Z+jets, and therefore they are combined in order 
to treat them as fully correlated. The normalisation 
of this sample is allowed to vary freely in the fit. 
Independent validation of the Z-|-jets background is 
carried out in two-lepton control regions. The agree- 
ment is good between data and MC simulation in 
both normalisation and shape. 

— Diboson WW. This sample is constrained with a 
30 % cross-section systematic uncertainty. The con- 
straint is necessary because of the relatively small 
contribution of the sample in most signal and con- 
trol regions and because no W^W^-dominated control 
region can be constructed; if the background were 
allowed to vary freely, then the fit may find a mini- 
mum with an unreasonably large or small contribu- 
tion from diboson WW events and hide some other 
effect with an artificial WW normalisation. 

— Charge flip. Charge mis-identification can occur due 
to physical effects, like lepton Bremsstrahlung, and 
detector effects, especially for high-px leptons with 
almost straight tracks. These effects generate back- 
ground in the same-sign dielectron and electron- 
muon channels. This background is negligable in the 
dimuon channel, where the contribution from both 
physical and detector effects is far smaller. The elec- 
tron charge-flip rate is measured as a function of rj in 



the data 68 , allowing MC simulation to model the 



lesser dependence on px- These charge-flip rates are 
applied to opposite-sign MC simulation events, pro- 
viding an estimate of the overall contribution from 
charge flip in these channels. The electron is ad- 
ditionally shifted and smeared to mimic the effect 
of charge mis-identification. This shift in the px is 
propagated through to the razor variables. The un- 



certainty from the charge flip probabilities domi- 
nates the uncertainty of this background. 
Fake leptons. The multijet background in the one- 
lepton signal regions, as well as the H^-|-jets, semi- 
leptonic tt, and multijet background in the two- 
lepton signal regions, comes predominantly from had- 
rons faking electrons and muons. This background is 
estimated using the "matrix method" 13,68|, using 
the number of baseline leptons not passing signal 
lepton requirements. The efficiency for a real lepton 
passing the baseline lepton requirements to pass the 
signal lepton requirements is estimated using Z MC 
simulation events. The rejection rate for fake leptons 
is estimated in data, using samples enriched in fake 
leptons. For electrons, the factors are derived and 
applied separately for inclusive samples of events 
and samples requiring a 6-tagged jet. Because this 
background accounts for all fake background, MC 
events in the one-lepton (two-lepton) channels are 
required to have at least one (two) prompt lepton(s) 
from a r lepton, W boson, Z boson, or sparticle. The 
uncertainty on this background estimate has a sta- 
tistical component from the number of events in the 
control region and a systematic component from the 
uncertainty on the scale factors. 
Some fraction of the events with same-sign, baseline 
leptons in the data may be due to charge flip. Thus, 
the matrix method overestimates somewhat the fake 
lepton background in the dilepton channels. In order 
to correct for this overlap, opposite-sign events in 
data containing baseline leptons that do not pass the 
signal lepton requirements are used. Each event is 
assigned a weight representing the likelihood of that 
event being subject to charge mis-identification. The 
weighted events are then presented as a negative 
component to the same-sign fake background distri- 
bution, such that the contribution to the same sign 
fake background from originally oppositely charged 
leptons is subtracted. 

Multijets in zero-lepton channel. Two specific con- 
trol regions constrain this background, and its nor- 
malisation is allowed to vary freely in the fit. Several 
different approaches are used to cross-check this es- 
timate. Prescaled single jet triggers are used to con- 
struct independent multijet-enriched control regions 
at low M'j^ that is free from the inefficiency of the 
i?x-based trigger. The observed number of events in 
this region are then projected into the signal region 
using transfer factors from MC simulation. Alter- 
natively, in order to model the mis-measurement of 
jets in the calorimeter, jets in events collected with 
these single-jet triggers are smeared according to re- 
sponse functions estimated using data |10J . Both of 
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Table 2 Background control, signal, and validation regions. All signal regions include the overflow in the highest bin. "N/A" 
means that there is no requirement. The binning of the validation regions does not affect the results, so they are listed as 
"N/A." 

Name Leptons 6-jets W,jets R Range Range Number of bins 



Control Regions 



Had. 6-veto Multijet 


leptons 


= 


> 5 


0.3 < i? < 0.4 


800 < Af' 


< 2000 GeV 


12 in A/^ 


Had. 6-tag Multijet 


leptons 


> 


> 5 


0.2 < i? < 0.3 


1000 <M'< 2000 GeV 


10 in Af^ 


e W+jets 


1 electron 


= 


> 5 


< i? < 0.7 


300 < M'j^ 


< 400 GeV 


7 in ii 


A* VK+jets 


1 muon 


= 


> 5 


< i? < 0.7 


300 < M'^ 


< 400 GeV 


7 in ii 


e tt 


1 electron 


> 


> 5 


< i? < 0.7 


400 < M^j 


< 650 GeV 


7 in J? 


fi tt 


1 muon 


> 


> 5 


< i? < 0.7 


400 < M'^ 


< 650 GeV 


7 in i? 


ee tt 


2 OS electrons 


> 


N/A 


0.2 < i? < 0.3 


400 < M'^ 


< 1000 GeV 


6 in M'j^ 


fiji ti 


2 OS muons 


> 


N/A 


0.2 < _R < 0.3 


400 < M'^ 


< 1000 GeV 


6 in M'j^ 


efi tt 


2 OS OF leptons 


> 


N/A 


R < 0.3 


400 < M'j^ 


< 1200 GeV 


8 in M'j^ 


ee Z 


2 OS electrons 


N/A 


N/A 


R < 0.4 


600 < M'j^ 


< 1500 GeV 


9 in M'j^ 


fifi Z 


2 OS muons 


N/A 


N/A 


R < 0.4 


600 < M'j^ 


< 1500 GeV 


9 in M'j^ 


ee Charge flip 


2 SS electrons 


N/A 


N/A 


R < 0.25 


500 < M^j 


< 1200 GeV 


7 in M'j^ 


Signal Regions 


Had. 6- veto 


leptons 


= 


> 5 


R > 0.70 


600 < M^j 


< 1200 GeV 


3 in M'j^ 


Had. b-tag 


leptons 


> 


> 5 


R > 0.40 


900 < M'^ 


< 1500 GeV 


3 in M'j^ 


e fa-veto 


1 electron 


= 


> 5 


R > 0.55 


500 < M' 


< 1000 GeV 


3 in M'j^ 


e 6-tag 


1 electron 


> 


> 5 


R > 0.35 


1000 < M' < 1600 GeV 


6 in M'j^ 


11 b-veto 


1 muon 


= 


> 5 


R > 0.55 


500 < M' 


< 1000 GeV 


3 in M'j^ 


fi b-tag 


1 muon 


> 


> 5 


R > 0.35 


1000 < M' < 1400 GeV 


4 in M'j^ 


OS-ee 


2 OS electrons 


N/A 


N/A 


R > 0.40 


600 < M'^ 


< 1000 GeV 


4 in M'j^ 


OS-^J,^l 


2 OS muons 


N/A 


N/A 


R > 0.40 


600 < M'^ 


< 1000 GeV 


4 in M'j^ 


SS-ee 


2 SS electrons 


N/A 


N/A 


R > 0.25 


500 < M'^ 


< 900 GeV 


4 in M'j^ 


SS-fifj, 


2 SS muons 


N/A 


N/A 


R > 0.25 


500 < M'j^ 


< 900 GeV 


4 in M'j^ 


OS-e/i 


2 OS OF leptons 


N/A 


N/A 


R > 0.40 


600 < M^j 


< 1000 GeV 


4 in M'j^ 


SS-e/x 


2 SS OF leptons 


N/A 


N/A 


R > 0.25 


500 < M'^ 


< 900 GeV 


4 in Af^ 


Validation Regions 


Had. b-veto Multijet 


leptons 


= 


> 5 


0.4 < i? < 0.6 


800 < M' 


< 2000 GeV 


N/A 


Had. b-tag Multijet 


leptons 


> 


> 5 


0.3 < i? < 0.4 


1100 < M' < 2000 GeV 


N/A 


1-lep b-veto W+jets 


1 lepton 


= 


> 5 


N/A 


400 < M^j 


< 550 GeV 


N/A 


1-lep b-tag tt 


1 lepton 


> 


> 5 


N/A 


700 < M' 


< 850 GeV 


N/A 


OS-ee/ fi^ tt 


2 OS SF leptons 


iO 


N/A 


0.3 < i? < 0.4 


400 GeV < M' 


N/A 


OS-efi ti 


2 OS OF leptons 


iO 


N/A 


0.3 < i? < 0.4 


N/A 




N/A 



these methods result in an estimate consistent with 
that derived in the main fit. 

The systematic effects included as nuisance param- 
eters in the fit are: the jet energy scale and resolu- 
tion uncertainties; 6-tagging uncertainties; uncertainty 
on the MC simulation modelling of the JVF; the ad- 
ditional cross-section uncertainty on the production of 
heavy flavour in association with a vector boson; the 
uncertainties on trigger efficiency and matching and re- 
construction efficiency; a systematic uncertainty on the 
re- weighting of the VF-boson p^; uncertainties on the 
missing transverse momentum pile-up dependence and 
the calibration of energy not associated with an object 
in the event; the matrix method statistical and sys- 
tematic uncertainties; the charge flip systematic uncer- 
tainties; the diboson WW shape systematic uncertainty 
taken from comparing Herwig to Alpgen. Where the 
systematic uncertainties affect object definitions, cor- 
rections are propagated to the missing transverse mo- 



mentum and razor variable calculations. The effects of 
other uncertainties on the final results are negligible. 
These uncertainties affect the signal yield and shape 
in the signal regions, as well as the allowed variation in 
signal-region background estimates after the control re- 
gion constraints. In most signal regions, the jet energy 
scale uncertainty is the dominant experimental uncer- 
tainty (from 10% to 25%). 



7 Background fit 

Figure [2] shows the distributions of Af^ and jet multi- 
plicity in the zero-lepton multijet control region with 
a 6-tagged jet requirement, with results from the fit 
to the control regions overlayed. By design, the multi- 
jet background is dominant in these regions. The small 
contribution from tt and W-|-jets backgrounds are con- 
strained by other control regions in the simultaneous 
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Fig. 1 A visual representation of the zero-lepton (top), one-lepton (middle) signal regions, and two-lepton (bottom) control 
(CR), validation (VR), and signal (SR) regions. 



fit. The hatched area indicates the total systematic un- 
certainty after the constraints imposed by the fit. 

The distributions of R and jet multiplicty for the 
backgrounds after the control region fit in the W 
/iz^+jets control region are sho-wn in Fig.|3] The control 
region at lo-w R is dominated by fake backgrounds, and 
at moderate-to-high R they are dominated by VF-|-jets. 



The use of an alternate control region "with a cut on 
transverse mass, -which significantly reduces the fake 
contribution, results in a negligible change in the final 
search results. 

Figure |4] and [5] sho"w the one-lepton and t-wo-lepton 
ti control regions, respectively. The fit reduces the nor- 
malisation of the tt background in the one-lepton con- 
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Fig. 2 The distribution of (left) and the number of jets with pT > 30 GeV (right) in the multijet control region with 
a 6-tagged jet requirement (dots with error bars), the expectation from the control region fit for various backgrounds (filled), 
and the systematic uncertainty (hatched). 




R N jets (p_|_>30 GeV) 

Fig. 3 The distribution of R (left) and the number of jets with pT > 30 GeV (right) in the W — )■ ^u+jets control region (dots 
with error bars), the expectation from the control region fit for various backgrounds (filled), and the systematic uncertainty 
(hatched). 



trol region by approximately 15-20% with respect to 
the unmodified expectation from MC simulation. Tliis 
shift predominantly affects the semi-leptonic tt back- 
ground. In the two-lepton analysis, there is a signifi- 
cant contribution to the background expectation from 
Z-boson events with heavy flavour, particularly at low 
^miss^ The lowest M}j bin shows the most significant 
disagreement, which demonstrates the importance of 
shape profiling by binning the control regions. Although 
in that lowest bin, particularly in the two-muon chan- 
nel, the MC simulation underestimates the amount of 
data, a single-binned normalisation of the background 
would result in an overestimation of the background at 
high M'j^. The distributions of missing transverse mo- 



mentum are also shown for the two-lepton control re- 
gions. 

The distributions of Af^ for the two Z+jets control 
regions are shown in Fig. [6j In both the electron and 
muon channel, good agreement after the control region 
fit is visible. 

The charge flip background is significant in the same- 
sign two-electron channel. Figure |6] also shows the dis- 
tribution of M'j^ in a charge flip enriched control region. 
There remain significant uncertainties on the background 
even after the control region fit, since it is dominated by 
charge flip and fake leptons, both of which have large 
systematic uncertainties associated with them. The dis- 
tribution of dilepton mass is also shown. 
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Fig. 4 The distribution of R in the one-lepton tt control regions (dots with error bars), the expectation after the control 
region fit for various backgrounds (filled), and the systematic uncertainty (hatched). 



Table 3 Results of the background-only fit to the control regions in the zero- and one-lepton control regions, for an integrated 
luminosity of 4.7 fb~^. The errors shown are the statistical plus systematic uncertainties. 



Control region 


Had. 6-veto Multijet 


Had. 6-tag Multijet 


e VF+jets 


fi VF+jets 


e ti 


^ tt 


Observed events 


1032 


2153 


1833 


1413 


3783 


3479 


Fitted background events 


1030 ± 30 


2150 ± 50 


1840 ± 40 


1410 ± 30 


3820 ± 60 


3470 ± 50 


Fitted top events 


21 ±7 


170 ± 19 


280 ± 30 


290 ± 30 


2800 ± 60 


2800 ± 60 


Fitted W/Z events 


90± 10 


26 ±4 


670 ± 40 


690 ± 50 


210 ± 20 


240 ± 30 


Fitted WW diboson events 


0.54 ±0.18 


0.14 ±0.05 


4.2 ± 1.8 


4.5 ± 1.9 


1.2 ±0.5 


1.0 ±0.4 


Fitted multijet events 


920 ± 30 


1960 ± 50 


0±0 


0±0 


0±0 


0±0 


Fitted charge flip events 


0±0 


0±0 


0±0 


0±0 


0±0 


0±0 


Fitted fake lepton events 


0±0 


0±0 


890 ± 50 


430 ± 40 


810 ± 70 


440 ± 60 



Expected events 


990 


2670 


2110 


1560 


4300 


3790 


MC exp. top events 


52 


245 


450 


470 


3300 


3250 


MC exp. W/Z events 


110 


28 


740 


760 


200 


220 


MC exp. WW diboson events 


0.61 


0.18 


4.6 


4.6 


1.4 


1.1 


MC exp. multijet events 


830 


2400 














Charge flip events (estimated from data) 




















Fake lepton events (estimated from data) 








910 


330 


800 


310 



The contributions to each of the control regions be- 
fore and after the fit to the control regions are shown 
in Tables [3] and m 

Various tests of the fit are carried out in order to 
ensure its stability. As a test of the multijet background 
constraint and the validity of fitting the M'j^ distribu- 
tions in those control regions, the control region fit is 
instead performed in the number of jets with px > 
30 GeV. The pt cut is raised from the baseline selec- 
tion to make the fit less sensitive to pile-up effects. The 
expectation for the multijet background in the signal 
regions is consistent with the main result. 

The yields and shapes of the validation regions show 
good agreement with the Standard Model expectation. 
The significance of the deviation of the observation from 



the expectation in each of the signal and validation re- 
gions are shown in Figure [7] There is some tension be- 
tween the same-flavour and opposite-flavour dilepton tt 
validation regions, but there is no indication of a sys- 
tematic mis-modelling of any of the major backgrounds. 
The yields of all validation regions are within 1.2a of 
the SM expectations. 

The numbers of expected events in each signal re- 
gion before and after the fit to the control regions are 
shown in Tables [5] and [6j Additionally, the probability 
(po-value) that a background-only pseudo-experiment 
is more signal-like than observed is given for each in- 
dividual signal region. To obtain these po-values, the 
fit in the signal region proceeds in the same way as the 
control-region-only fit, except that the number of events 
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Fig. 5 The distribution of M'^ (left) and missing transverse momentum (right) in the two-lepton tt control regions (dots 
with error bars), the expectation after the control region fit for various backgrounds (filled), and the systematic uncertainty 
(hatched). 
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Table 4 Results of the background-only fit to the control regions in the two-lepton control regions, for an integrated luminosity 
of 4.7 fb~^. The errors shown are the statistical plus systematic uncertainties. 



Control region 


ee tt 


tt 


ep, tt 


ee Z 


/i/i Z 


ee Charge flip 


Observed events 


272 


347 


1340 


3688 


4579 


183 


Fitted background events 


277 ± 14 


310 ± 10. 


1320 ± 30 


3670 ± 60 


4590 ± 70 


183 ± 13 


Fitted top events 
Fitted W jZ events 
Fitted WW diboson events 
Fitted multijet events 
Fitted charge flip events 
Fitted fake lepton events 


198 ±7 
45 ±4 
0.22 ±0.08 
0±0 
0±0 
34± 15 


237 ±8 
51 ± 5 
0.10 ± 0.15 
0±0 
0±0 
22 ±8 


1090 ± 30 
3.5 ± 0.3 
1.3 ± 0.5 
0±0 
0±0 
220 ± 40 


220 ±9 
3090 ± 90 
6± 3 
0±0 
0±0 
360 ± 100 


281 ± 11 
4220 ± 80 
8± 5 
0±0 
0±0 
80 ±50 


0.104 ±0.011 
1.06 ± 0.11 
1.2 ±0.6 
0± 
94 ± 14 
87± 19 



Expected events 


305 


336 


1340 


3920 


5050 


148 


MC exp. top events 


225 


276 


1220 


278 


357 


0.094 


MC exp. W jZ events 


41 


47 


3.1 


3360 


4600 


1.14 


MC exp. WW diboson events 


0.21 


0.09 


1.2 


6 


8 


1.2 


MC exp. multijet events 




















Charge flip events (estimated from data) 

















94 


Fake lepton events (estimated from data) 


39 


13 


120 


270 


80 


51 



Table 5 Results of the background-only fit to the control regions in the zero- and one-lepton signal regions, for an integrated 
luminosity of 4.7 fb~^. The errors shown are the statistical plus systematic uncertainties. The po-values and significances are 
given for single-bin signal regions with somewhat tighter cuts, along with the 95 % Confidence Level upper limit on the 
number of and cross-section for non-Standard-Model production within each signal region. 



Signal region 


Had. b-veto 


Had. 6-tag 


e fe-veto 


e b-tag 


6-veto 


{i 6-tag 


Observed events 


4 


30 


6 


13 


9 


4 


Fitted background events 


5.5 ± 1.5 


39 ±7 


10 ±2 


6.6 ± 1.7 


5.5 ± 1.7 


4.4 ± 1.3 


Fitted top events 


0.40 ±0.14 


21 ±3 


2.7 ±0.9 


5.0 ± 1.3 


1.7 ±0.6 


3.7± 1.1 


Fitted W jZ events 


4.9 ± 1.3 


3.8 ±0.7 


7.2 ± 1.7 


1.2 ±0.5 


3.8 ± 1.3 


0.6 ±0.5 


Fitted WW diboson events 


0.03 ± 0.02 


0.029 ±0.010 


0.01 ± 0.02 


0.000 ± 0.009 


0.001 ± 0.008 


0.010 ± 0.005 


Fitted multijet events 


0.25 ±0.10 


14 ±5 


0±0 


0±0 


0±0 


0±0 


Fitted charge flip events 


0±0 


0±0 


0±0 


0±0 


0±0 


0±0 


Fitted fake lepton events 


0±0 


0±0 


0.3 ±0.3 


0.5 ±0.6 


0±0 


0±0 




Expected events 


6.7 


55 


14 


8.5 


9.5 


5.1 


MC exp. top events 


0.88 


30 


5.7 


6.3 


3.4 


4.6 


MC exp. W jZ events 


5.6 


4.0 


8.5 


1.8 


6.1 


0.5 


MC exp. WW diboson events 


0.04 


0.046 


0.01 


0.000 


0.012 


0.010 


MC exp. multijet events 


0.20 


21 














Charge flip events (estimated from data) 




















Fake lepton events (estimated from data) 








0.3 


0.5 










Tight M^j cut (GeV) 


600 


1100 


600 


1100 


600 


1100 


Observed events 


4 


5 


5 


6 


2 


4 


Background events 


6.2 ± 1.8 


13 ±3 


5.3 ± 1.6 


2.4 ± 1.0 


2.4 ± 1.0 


1.9 ±0.8 



po-value (Gauss, cr) 0.72 (-0.57) 0.91 (-1.35) 0.53 (-0.07) 0.07 (1.50) 0.54 (-0.10) 0.16 (0.98) 

Upper limit on A^bsm 5.2 (6.3|^'^) 6.5 (9.3|f,^) 6.3 (6.4|!^^^) 9.0 (5.5||'^) 4.4 (4.5|^-J) 6.8 (4.8|J'^) 

Upper limit on a (fb) 1-1 (l^ag) 1-4 (S-Oj?;^) 1.3 (1.4|g:°) 1.9 (1.2|J:|) 0.9 (l.ojj;^) 1.4 (l.ojj;?) 
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Fig. 6 Top, the distribution of M'^ in the Z+jets control regions. Bottom, the distribution of (left) and dilepton mass 
(right) in the charge flip control region (dots with error bars), the expectation after the control region fit for various backgrounds 
(filled), and the systematic uncertainty (hatched). 



observed in the signal region is included as an input 
to the fit. Then, an additional parameter for the non- 
Standard-Model signal strength, constrained to be non- 
negative, is fitted. The shape of the distributions in the 
signal region is neglected in this fit. Therefore, in order 
to provide tighter constraints on non-Standard-Model 
production, in some of the high-count signal regions 
the M'j^ requirements are tightened. The significance of 
the excess is given, along with the model-independent 
upper limit on the number of events and cross-section 
times acceptance times efficiency from non-Standard- 
Model production. The distributions in all signal re- 
gions as a function of Af^ of background expectations 
after the fit to the control region has been performed are 
shown in Figs. [8] and |9] No significant deviations from 
the expected background are found. The most signifi- 
cant excess is 1.50 standard deviations from the expec- 
tation, in the one electron, 5-tagged jet signal region. 



8 Exclusion results 

Using these signal regions, the CLs [69j prescription is 
applied to find 95% Confidence Level (CL) one-sided 
limits on the production of SUSY events in various 
models. A likelihood is constructed, taking into account 
signal shape information provided by the binning of the 
signal regions. All fitted nuisance parameters, with their 
correlations, are included in the likelihood. Because the 
typical AI'j^ of the signal may vary across a signal grid, 
the use of shape information results in an observed ex- 
clusion that is not consistently above or below the ex- 
pected. The observed limits for the separate zero-, one-, 
and two-lepton signal regions are also constructed addi- 
tionally, with all control regions included as constraints. 



Figure 10 shows exclusion contours for a simplified 
model with gluino pair-production, where the gluinos 
decay to a chargino and two quarks and the chargino 
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Table 6 Results of the background-only fit to the control regions in the two-lepton signal regions, for an integrated luminosity 
of 4.7 fb~^. The errors shown are the statistical plus systematic uncertainties. The po-values and significances are given for 
single-bin signal regions, along with the 95 % Confidence Level upper limit on the number of and cross-section for non-Standard- 
Model production within each signal region. 



Signal region 


OS-ee 




SS-ee 




OS-e/i 


SS-e/i 


Observed events 


10 


15 


11 


8 


18 


18 


Fitted background events 


12 ± 2 


13 ±2 


6±4 


4± 3 


20. ±3 


14 ± 8 


Fitted top events 


10.2 ± 1.5 


10.7 ± 1.6 


0.12 ±0.04 


0.39 ±0.17 


19 ±2 


0.7 ± 0.2 


Fitted W jZ events 


0.54 ±0.10 


0.6 ±0.2 


0.16 ±0.04 


0.10 ±0.04 


0.26 ± 0.04 


0.33 ±0.07 


Fitted WW diboson events 


0.4 ±0.4 


0.4 ±0.5 


0.6 ± 0.3 


0.6 ±0.5 


0.5 ± 1.0 


1.2 ±0.7 


Fitted multijct events 


0±0 


0±0 


0± 


0± 


0±0 


0±0 


Fitted charge flip events 


0±0 


0±0 


1.6 ± 0.4 


0± 


0±0 


1.1 ±0.2 


Fitted fake Icpton events 


1.2 ± 1.3 


1.3 ± 1.1 


3±4 


3± 3 


0.6 ±0.6 


10. ± 8 




Expected events 


15 


16 


6 


5 


24 


14 


MC exp. top events 


13.1 


14.7 


0.13 


0.49 


23 


0.6 


MC exp. W jZ events 


0.67 


0.4 


0.19 


0.21 


0.27 


0.36 


MC exp. WW diboson events 


0.4 


0.4 


0.7 


0.7 


0.5 


1.2 


MC exp. multijet events 




















Charge flip events (estimated from data) 








1.6 








1.0 


Fake lepton events (estimated from data) 


1.2 


1.3 


3 


3 


0.6 


11 




po-value (Gauss, cr) 


0.71 (-0.56) 


0.32 (0.46) 


0.15 (1.05) 


0.18 (0.93) 


0.68 (-0.48) 


0.29 (0.54) 


Upper limit on A^bsm 


7.3 (8.8l^^/) 


11.1 (9.4j-/) 


14.0 (10.2t^V) 


11.4 {8.0j;V) 


9.4 (ll.lI^V) 


17.7 (i4.9i?;;:^) 


Upper limit on a (fb) 


1.6 (i.gt?-^) 


2.4 (2.0j?;») 


3.0 (2.2t-^;i) 


2.4 


2.0 (2.4j-^;^,) 


3.8 (3.2t*;*) 
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Fig. 7 Pull distributions of the numbers of events in the validation regions (VR) and signal regions (SR). The filled (dashed) 
bars show the agreement after (before) the background-only fit to the control regions has been performed. 
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Fig. 8 The all-hadronic (top) and one-lepton (bottom) signal regions with a 6-tagged jet veto (left) and requirement (right), 
after the fit to the control regions has been performed. 



subsequently decays to a W boson and the LSP. Two 
planes are shown for this simplified model. The first 
fixes the chargino mass to be exactly between the LSP 
and gluino mass and shows the exclusion in the gluino- 
mass-LSP-mass plane. The production cross-section falls 
smoothly and exponentially with mheavy, while M'jj and 
therefore the acceptance times efficiency for a signal re- 
gion typically rises with the mass splitting, mheavy ^ 
TTiLSP. In the second plane, the LSP mass is fixed to 
60 GeV and the exclusion is shown in the gluino mass-a; 
plane, where x = (mchargino - ^lsp)/ (m-giuino - "Ilsp)- 
The zero- and one-lepton signal regions with a fe-tagged 
jet requirement do not contribute to the exclusion be- 
cause these simplified models have only light quarks in 
the matrix element final state. 

At high X, although the leptons have high-px, the 
larger branching fraction of the W to quarks allows the 
zero-lepton channel to dominate. At moderate x, the 
leptons allow better discrimination between signal and 



background in the one- and two-lepton channels. At 
low X, the leptons have too low pt, and the zero-lepton 
channel again dominates. At high x, the limit set by 
this analysis exceeds somewhat that of the dedicated 



0-lepton and 1-lepton ATLAS searches 10 13 , while at 



low X the limit is weaker. This asymmetry is produced 
by the difference in kinematics in these two regions of 
the plane. At high x, the charginos are almost at rest in 
the lab frame, and the event topology is dominated by 
a two-body decay, Xi ~^ ^^'''Xi- At low x, on the other 
hand, the chargino is highly boosted, and the topology 
is dominated by a three-body decay, g — )> qq'x^ ■ Thus, 
the high-x events have a higher typical R than the low- 
X events, and the two have approximately the same 
distribution. 



Figure 11 shows exclusion contours in simplified mod- 
els with gluino pair-production, where the gluinos decay 
to the LSP via the emission of a tt pair. The exclu- 
sion is presented in the gluino mass-LSP mass plane. 
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Fig. 9 The two-lepton signal regions for same-flavour (top) and opposite-flavour (bottom) leptons of the same sign (left) and 
opposite sign (right), after the fit to the control regions has been performed. 



and, since all top quarks are required to be on-shell, 
only points with mgiuino > JTiLSP + 2 x mtop are con- 
sidered. The zero- and one-lepton signal regions with 
a 6-tagged jet veto do not contribute to the exclusion, 
because these models include four top quarks per event. 
At small mass splitting, the limits here are somewhat 
stronger than the ATLAS multijet analysis though 
both are exceeded by those of the dedicated multi-6- 
jet analysis p^. The three 5-tagged jet requirement 
suppresses the background substantially while preserv- 
ing the signal acceptance because of the four tops in 
the event. The combined limit on LSP mass falls more 
quickly than that of the multi-6-jet analysis because 
M'j^ is proportional to the mass splitting in the event, 
here mgiuino ~ m,LSP- The zero-lepton razor analysis is 
limited in this case by the use of the Ht trigger, which 
was chosen to avoid a bias in the distribution. 



/X > 0. At low mo, where squark pair-production is dom- 
inant, the zero-lepton channel dominates the exclusion, 
although it is affected somewhat by the jet multiplicity 
requirement that is not applied in the dedicated signal 



region of Ref. 10 , which therefore has more stringent 



Finally, Figure 12 shows exclusion contours in a 
plane of MSUGRA with tan(/3) = 10, = 0, and 



limits. The leptonic channels enter at high mo, par- 
ticularly where longer decay chains are common. The 
robustness of the individual limits have also been cross 
checked by removing some of the control regions. For 
example, removing the zero- and one-lepton control re- 
gions from the calculation of the two-lepton limit, the 
MSUGRA limit changes by less than 20 GeV in mi/2- 
These limits are consistent with those of earlier ATLAS 
analyses [l0][TT}[l3] in the mgiuino ~ msquark region. In 
this region, the single mass-splitting scale of the main 
strong production modes should produce a somewhat 
sharper peak in M'j^, allowing an improved limit in the 
shape fit. At large mo, the high M'j^ requirement of 
the all-hadronic signal regions, resulting from the Ht 
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Fig. 10 Exclusion in a simplified model with gluino pair-production, where the gluinos decay to a chargino via the emission 
of two quarks and the chargino decays to the LSP and a W boson. Top, for the chargino mass exactly between the gluino and 
LSP mass, in the gluino mass-LSP mass plane, and bottom, for the LSP mass fixed to 60 GeV, in the gluino mass-x plane, 
with X = (mchargino — "^LSP )/('Tigiuino — 'T^LSp)- The cxclusiou is shown for the combination, as well as for each individual 
channel. Overlaid is the exclusion result published in Ref. ,13,. 
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Fig. 11 Exclusion in a simplified model with gluino pair-production, where the gluinos decay to the LSP via the emission of 
a tt pair. The exclusion is shown for the combination, as well as for each individual channel. Overlaid is the exclusion result 
published in Ref. [12]. 
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Fig. 12 Exclusion in a plane of the constrained minimal supersymmetric model. The exclusion is shown for the combination, 
as well as for each individual channel. Overlaid is the exclusion result published in Ref. [10]. 
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trigger use, produce a somewhat weaker limit than the 



ATLAS multijet analysis 11 



The complementarity of a search using razor vari- 
ables can be quantified by studying the overlap of the 
signal regions with the dedicated searches. Various sig- 
nal models have been studied to understand this over- 
lap, including both simplified models and full SUSY 
production models. The overlap between the signal re- 
gions presented here and other searches in ATLAS 10 



11 13 Iff is typically 10-50%, with similar overlaps 



in the data. The signal regions of this search access 
kinematic regions that are different from those of the 
standard searches. In simplified models in particular, 
the overlap between the dominant signal regions in the 
standard ATLAS analyses and the signal regions pre- 
sented here is below 10-15 %. Thus, the regions of SUSY 
parameter space and kinematic phase space excluded by 
this search complement those excluded by earlier AT- 
LAS searches using the same data sample. 



9 Summary 

A search for supersymmetry including final states with 
zero, one, and two leptons, with and without fe-tagged 
jets, in 4.7 fb~^ of ^/s = 7 TeV pp collisions has been 
presented. Mutually exclusive signal regions exploiting 
these final states are combined with the use of variables 
that include both transverse and longitudinal event in- 
formation. No significant excess of events beyond the 
Standard Model background expectation was observed 
in any signal region. Fiducial cross section upper lim- 
its on the production of new physics beyond the Stan- 
dard Model are shown. Exclusion contours at 95 % CL 
are provided in SUSY-inspired simplified models and 
in the constrained minimal supersymmetric extension 
of the Standard Model. 
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M. Medinnis*^, S. Meehan^i, R. Meera-Lebbai^^S T. Meguro"^, S. Mehlhase^'^, A. Mehta^^, K. Meier^^a, 
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A. Mengarelli^O'^'^ob, g. Menke^^, E. Meonii^i, K.M. Mercurio", P. Mermod^'^, L. Merolai°2a,i02b^ q Meroni^S'*, 

F. S. Merritt^S H. Merritt^^^, A. Messina^O'^', J. Metcalfe^^, A.S. Mete^*^^ C. Meyer«\ C. Meyer^^, 
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R. Prabhu^'^, P. Pralavorio^', A. Pranko^^ S. Prasad'", R. Pravahan25, S. Prell^', K. Pretzl^'^, D. Price'^", 
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D. Silvcrstcini43, S.B. Silvcrstcin"6a^ Simaki^e, Q. Simardi'^^ Lj. Siniici-^'\ S. Simion"-\ E. Simioni**i, 

B. Simmons^^ R. SimonieUo^^'^'S^'', M. Simonyan^e, R Sinervol5^ N.B. Sinev""*, V. Sipica"\ G. Siragusa^^"^ 
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T. Sumida^^, X. Sun-"^-\ J.E. Sundcrmann*^ K. Suruliz^'^^, G. Susinno-^^-'^'^^b^ M.R. Suttoni^^, Y. Suzuki6'\ 
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